Superconductivity of calcium under the pressure at 120 GPa 
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The properties of the superconducting state in calcium under the pressure at 120 GPa were an- 
alyzed. By using the imaginary axis Eliashberg equations it has been shown, that the Coulomb 
pseudopotential reaches the high value equal to 0.215. In the considered case, the critical temper- 
ature is not properly described by the Allen-Dynes formula and it should be calculated with an 
use of the modified expression. In the paper the exact solutions of the Eliashberg equations on 
the real axis were also obtained. On this basis it was stated, that the effective potential of the 
electron-electron interaction is attractive for the frequencies lower or equal to the maximum phonon 
frequency. Then, the dimensionless parameter 2A (0) /fcsTc ~ 4.10 was calculated. In the last step 
it has been proven, that the ratio of the electron effective mass to the bare electron mass is high 
and reaches its maximum equal to 2.36 for the critical temperature. 

PACS numbers: 74.20.Fg, 74.25. Bt, 74.62.Fj 
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I. INTRODUCTION 

The first mention of the pressure-induced supercon- 
ducting state in calcium was given by Dunn and Bundy 
in 1981 f]\. In 1996 Okada et al. have obtained experi- 
mentally the dependence of the critical temperature (Tc) 
on the pressure (p) [2]. The authors stated, that the crit- 
ical temperature grows together with an increase of the 
pressure from the value of about 1 K for p = 50 GPa to 
the value ~ 20 K for 150 GPa. The results of Okada were 
verified in 2006 by Yabuuchi et al. The significant 
discrepancy between new and old data was observed (see 
Fig.[T]). Most probably the difference between the exper- 
imental results comes from the use of the samples of dis- 
tinct purity or from the different methods of the pressure 
measuring. The discussion of the considered problem the 
reader can find in the paper fsi]. 

According to Yabuuchi et al. the superconducting 
state in calcium can be observed in the following crys- 
tal structures: sc, Ca-IV and Ca-V. In the structure sc 
the critical temperature quickly increases together with 
the growth of the pressure (from the value of about 3 
K for p = 58 GPa to the value ~ 23 K for the pres- 
sure equal to 113 GPa). In the case of Ca-IV and Ca-V 
the increment of Tc is much slower. Although, with the 
further compression to 161 GPa the critical temperature 
increases to the value equal to 25 K (the highest observed 
value of the critical temperature for the simple metals) . 

The issue of the Tc value for the given pressure seems 
to be solved by Yabuuchi et al.. Unfortunately, the deter- 
mination of the proper sequence of the structural tran- 
sitions is still the tangled problem. On the one hand it 
is connected with the objective experimental difficulties, 
on the other hand the ab initio calculations are also bur- 
dened with the large uncertainty due to the small differ- 
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ence between the values of the enthalpy of the considered 
crystal structures [l]. 

Basing on the results presented in Q-@ it can be 
stated, that in the normal conditions calcium crystallizes 
in fee structure (Fig.[T] (A)). Under the influence of the 
pressure 20 GPa the structural transition fee-bee was ob- 
served; above 32 GPa calcium crystallizes in se structure. 
The crystal structure se is stable up to the pressure of 
about 113 GPa, where another structural change (Ca-IV) 
was discovered. The last transition to Ca-V structure was 
observed for the pressure of 139 GPa. It is worth men- 
tioning, that the crystal structures of the draft symbols 
Ca-IV and Ca-V were experimentally classified by Fuji- 
hisa et al. in the paper [10[. In particular, the structure 
Ca-IV is identified with P41212 and Ca-V with Cmea. 

Let us point the reader's attention toward the fact, 
that some recently published papers partially question 
above scheme of the structural transitions (see Fig. [1](B)- 
(D)) 0, 0. On the basis of the presented data it can 
be easily noticed, that the qualitative agreement between 
the obtained results exists only for the area of the low 
pressure's values (phases fee and bee). Above ~ 36 GPa 
the differences become significant. Thus at the present 
stage of the research, it is very difficult to specify which 
of the presented schemes is proper. 

From the physical point of view, the superconducting 
phase in calcium should have the most interesting prop- 
erties at the area of high pressures, where Tc > 20 K. 
For that reason, the thermodynamic parameters of su- 
perconducting state under the pressure at 120 GPa were 
determined. In this case Yin et al. suggest, that struc- 
tures: -P43212J Cmea and Pnma are showing almost iden- 
tical values of the enthalpy (Fig.[T] (C)) [8]. The paper 
of Ishikawa et al. predicts the appearance of the crystal 
structure Pnma (Fig.[T](D)) 
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FIG. 1: (A) The dependence of the critical temperature on 
the pressure for calcium. Experimental data: stars - S. Okada 
et al. [31, circles - T. Yabuuchi et al. [3]. Additionally, we 
present a sequence of the successive structural transitions pre- 
dicted by: (A) - T. Yabuuchi et al. 0] (see also the papers 
H-Sl), (B) - S. Arapan et al. 01, (C) - Z.P. Yin et al. d, 
(D) - T. Ishikawa et al. 



II. THE ELIASHBERG FORMALISM 

In the paper, the thermodynamic parameters will be 
calculated in the strict way in the framework of the 
Eliashberg formalism [Tl| . In particular, the critical 
value of the Coulomb pseudopotential {^q) and the pa- 
rameters in the modified Allen-Dynes formula are going 
to be determined with an use of the Eliashberg equa- 
tions defined on the imaginary axis [l3|-[i3]- The remain 
thermodynamic quantities (the order parameter and the 
electron effective mass) will be determined by using the 
Eliashberg equations in the mixed representation: the 
equations defined on both imaginary and real axis [Tsj . 



A. The Eliashberg equations on the imaginary axis 
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FIG. 2: (A) The order parameter on the imaginary axis for 
the selected values of the Coulomb pseudopotential (the first 
100 values of Am is presented). (B) The full dependence of 
the maximum value of the order parameter on the Coulomb 
pseudopotential. 
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where the symbol = (ia;„) denotes the order pa- 
rameter function and Zn = Z (iaj„) is the wave func- 
tion renormalization factor; n-th Matsubara frequency 
is given by the expression: aj„ = (tt//?) {2n — I), where 
P = {kBT)~^ and ks is the Boltzmann constant. The 
order parameter is defined as the ratio: A„ = 4>nlZn. 
The pairing kernel for the electron-phonon interaction is 
described with an use of the formula: 



A(z) = 2 



n 



(3) 



The Eliashberg function for calcium {a^F (il)) was deter- 
mined in the paper Q; the maximum phonon frequency 
(^^max) is equal to 61.68 meV. 

In the framework of the Eliashberg formalism, the de- 
pairing interaction between the electrons is taken into 
consideration with an aid of the Coulomb pseudopoten- 
tial fi* . The symbol 6 denotes the Heaviside unit function 
and uJc is the cut-off frequency; ujc — 3Qriia.x- 

The Eliashberg equations were solved for 2200 Matsub- 
ara frequencies (M = 1100) with an use of the iteration 
method, described in the papers fl^ and Tt]. In the 
considered case the solutions of the Eliashberg equations 
are stable for T > 5 K. 



The imaginary axis Eliashberg equations can be writ- 
ten in the following form [ll|: ^- The critical value of the Coulomb 



pseudopotential 

(j)^^ (1) The Coulomb pseudopotential is the second input pa- 

^ m=-M V + <^m ramcter in the Eliashberg equations. In general, its crit- 



TT ^ X{iuJn - ioJm) ~ {oJc - \i^m\) 
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FIG. 3: The value of the critical temperature as a function of 
the Coulomb pseudopotential. The filled circles represent the 
exact results obtained with an use of the Eliashberg equations. 
The solid line was achieved using the modified Allen-Dynes 
formula. The arrow shows the value — 0.215. The empty 
squares and triangles are related respectively to the classical 
Allen-Dynes expression and the McMillan formula. 



ical value may be theoretically determined from the first 
principles, however such types of calculations are very 
complicated and obtained results are usually burdened 
with the large error. According to the above, the param- 
eter /i^ will be selected in such way, that Tc determined 
on the basis of the Eliashberg equations equals the exper- 
imental value of the critical temperature. In the paper 
Tc = 24 K was assumed, in fair agreement with the re- 
sults obtained by Yabuuchi (see Fig.[T]). 

In order to succeed in above, the Eliashberg equations 
were solved in a strict way, for successive larger values 
of the Coulomb pseudopotential (T = Tc). Achieved 
results are presented in Fig. [5] (A), where the dependence 
of the order parameter on the Matsubara frequencies is 
being plotted. It can be easily noticed, that together 
with the increase of fi* the maximum value of the order 
parameter (Am=i) quickly decreases (Fig.[2](B)). On the 
basis of the plotted results, it has been stated, that /z^ = 
0.215. 

Let us notice, that the critical value of the Coulomb 
pseudopotential determined with an use of the Eliashberg 
equations significantly differs from the value /i^ calcu- 
lated with an aid of the classical formulas. In particular, 
the Allen-Dynes expression predicts: /z^ = 0.168 [18]; the 
estimations made on the basis of the McMillan formula 
are even more inaccurate {(1^—0.151) (Tgj . 



C. The critical temperature 

In the previous paragraph we have shown that the 
Allen-Dynes or McMillan formula can not be used in the 
case of calcium. For that reason, the values of the pa- 



rameters appearing in the classical Allen-Dynes formula 
were determined once again. In particular, we were bas- 
ing on the 200 exact values of Tc (/U*) determined with 
an use of the Eliashberg equations. The following result 
was received: 



ksTc = /i/2-p7^exp 
1.45 



-1.03(1 + A) 



\- fi* (1 -f G.06A) 



(4) 



where the symbol fi (/2) is the strong-coupling correction 
(shape correction) function: 
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The parameter uj2 denotes the second moment of the 
normalized weight function: uj2 = j /o^"""" dfla^F {fl) O. 
The quantity ujin is called the logarithmic phonon fre- 
quency and A is the electron-phonon coupling con- 



stant: Win 



exp 



2 f^i 

A Jo 



and A — 



2 /q^"''' dfl " ^^"^ . For calcium it was achieved: ^/uJ^ — 
29.81 nieV, win = 25.34 meV and A 1.3. In the paper 
we have changed also the parametrization of the func- 
tions Ai and A2. We have obtained: 



and 



Ao 



Ai = 2(l + 4.7Ai*) 



-0.085(1 - 150^i*) 



Win 
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In Fig. [3] we present a dependence of the critical tem- 
perature on the Coulomb pseudopotential. The results 
were obtained with an use of the Eliashberg equations, 
the modified and classical Allen-Dynes formula. Addi- 
tionally, the values of Tc calculated on the basis of the 
McMillan formula (/i = /2 = 1) are also plotted. It 
can be easily noticed, that the modified Allen-Dynes for- 
mula allows to very precisely reproduce the results of the 
Eliashberg theory in the whole considered range of /i*. 



D. The dependence of Am and Z„ 
temperature 



on the 



In the paragraph there are presented solutions of the 
Eliashberg equations on the imaginary axis for the tem- 
perature range from 5 K to Tc- In particular, in Fig. 3] 
(A) we have plotted the dependence of the order param- 
eter's values on the successive Matsubara frequencies. It 
can be easily noticed, that with the growth of the tem- 
perature the maximum of the order parameter's function 
{Am=i) is decreasing. Additionally, the half- width of the 
function becomes successively smaller. The last property 
means, that together with the temperature's growth less 
successive Matsubara frequencies commit a relevant con- 
tribution to the Eliashberg equations. Let us notice, that 
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FIG. 4: (A) The order parameter on the imaginary axis for 
the selected temperatures; first 100 values of Am is pre- 
sented. (B) The dependence of the maximum value of the 
order parameter on the temperature. The function can be 
parameterized with an use of the expression: Am^i {T) — 



Am^i (0) 
13 = 3.7. 



l-(^^y, where A™=i (0) = 4.16 meV and 



the dependence of the order parameter on temperature 
most conveniently can be traced after plotting function 
Am=i (T) (see Fig.|4] (B)). The obtained result proves, 
that the maximum value of the order parameter becomes 
saturated for the temperature equal to 5 K. For that rea- 
son the order parameter for T = 5 K can be used in the 
calculations of the energy gap at the temperature of zero 
Kelvin. 

The second solution of the Eliashberg equations has 
been shown in Fig.[S] (A). Also in this case, for the in- 
creasing values of the Matsubara frequencies the values of 
Zm are decreasing. However, in comparison with the or- 
der parameter, the decrease is much slower. The depen- 
dence of the maximum value of the wave function renor- 
malization factor on the temperature looks completely 
different than for Am=i (Fig. [5] (B)). It can be seen, that 
together with the growth of temperature Zrn=i has the 
shallow minimum and then slightly increases. 
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FIG. 5: (A) The wave function renormalization factor on the 
imaginary ax;is for the selected temperatures; first 200 values 
of Zm is presented. (B) The dependence of the maximum 
value of the wave function renormalization factor on the tem- 
perature. 



E. The Eliashberg equations in the mixed 
representation 



In order to exactly calculate the value of the energy 
gap at the temperature of zero Kelvin and electron ef- 
fective mass one should have the solutions of Eliashberg 
equations on the real axis. From the mathematical point 
of view the determination of such type function is a very 
complicated issue. Usually to do that, one uses the so- 
lutions on the imaginary axis and then analytically con- 
tinues them on the real axis |20j- However, the described 
procedure has a serious drawback e.g. the obtained func- 
tions are stable only for the low values of the frequencies 
(much lower than the maximum phonon frequency). For 
that reason, we have used a method of calculation which 
is deprived of the weak point mentioned above [l^. In 
particular, the approach is based on the transformation 
of the Eliashberg equations to the mixed representation, 
in which the order parameter function and the wave func- 
tion renormalization factor are defined on both real and 
imaginary axis. In the considered case the Eliashberg set 
takes the following form: 



5 



M 

(l){uj + iS) = - ^ [A - ioJm) - i^*d {<^c - l^^' 



and 

Z{uj + iS) 



in I du! a^F ( u! 



m=-M 
'0 



+ 00 

in I du! a^F 



N {oj ] + f {u} -oj 



(8) 



M 

^ rn = -M 



tn 



— / a^F { uj 

^ Jo 



r+oc 

ml , ' 2 7^ 

— / dui a F [UJ 



— cj + i5^ 



\J {lo — uj' Z"^ (a; — + i6) — (jp' {uj ~ lu' + iS) 
(f) (^uj + LU + iS 



\J{ijj + ijj' f Z2 (w + cj' + i6) - (j? {uj + uj' + iS) 



(9) 



^w — cj^Z^cj — cj + i6^ 



(w — cj' ) (oj — w' + iS) — 02 (oj — w' + i(5) 
cj + cj ] Z I OJ + OJ + iS 



^{uj + uj'f Z^ {uj + oj' + iS) ~ 41'^ {uj + oj' + iS) 








^^^^ 




^.^^^.^^ 




-•-Re[A(«))J 

^:^lm[A(m)l 




T=17.4 K 







10 20 30 40 50 60 

CO meV 



10 20 30 40 50 60 

CO meV 



12 






9 
6 






3 






^0 


-•- Re[4(«))] 
lm[A(io)] 

10*a¥(M) 




-3 


T=20.9 K 



10 20 30,40 50 60 
CO meV 




10 20 30,40 50 60 
CO meV 



FIG. 6: The dependence of the real and imaginary part of 
the order parameter on the frequency for the selected temper- 
atures. The rescaled Eliashberg function is also plotted. 



where symbols N [lS) and / (w) denote the statistical 
functions of bosons and fermions respectively. 



F. The order parameter on the real axis 

The Eliashberg equations in the mixed representation 
were solved for identical temperatures like the equations 
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FIG. 7: The values of A (lo) on the complex plane for the 
selected temperatures. 



on the imaginary axis. In Fig.|6] we have presented 
Re[A (w)] and Im[A (w)] for the range of the frequencies 
from to ilmax and selected temperatures. Additionally, 
the form of the Eliashberg function is also plotted. It can 
be easily noticed, that for the low frequencies the non- 
zero is only the real part of the order parameter. From 
the physical point of view it means the lack of the damp- 
ing effects, which are described by the imaginary part 
of the order parameter pTj . For higher frequencies the 
functions Re[A (lS)\ and Im[A (w)] are characterized by 
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the comphcated course which is very plainly correlated 
with the shape of the Eliashberg function. 

The full dependence of the order parameter on the tem- 
perature for the range of the frequencies from to Uc 
can be traced in the most convenient way on the com- 
plex plane (see Fig. [71). On the basis of the presented 
results it can be seen, that the values of the order pa- 
rameter are forming the characteristic spirals with radius 
that decreases together with the temperature growth. In 
Fig. [7] the lines with open symbols represent the solu- 
tions for u) G (0,r2niax); whereas solid lines correspond 
to the solutions for w S (fimax, "^c)- In the first case the 
condition Re[A (w)] > is fulfilled with a good approx- 
imation and in the second case Re[A (w)] < occurs. 
Let us notice, that the real part of the order parameter 
is connected with the effecti ve p airing potential for the 
electron-electron interaction [21|. Thus, the obtained re- 
sult means, that only in the range of the frequencies for 
which the Eliashberg function is defined, the effective po- 
tential for the electron-electron interaction is attractive. 

Below on the basis of the solution for T = 5 K we have 
calculated the value of the order parameter at the tem- 
perature of zero Kelvin (A (0)). In particular to achieve 
that, the following equation was used: 



A (T) = Re [A = A (T))] 



(10) 



As a result it was obtained: A(0) = 4.24 meV. The famil- 
iarity with the energy gap (2A(0)) allows to determine 
the dimensionless ratio i?i = , which for the calcium 
is equal to 4.10. Let us notice, that in the framework 
of the BCS theory the parameter Ri takes considerably 
smaller value that is equal to 3.53 [2^ . 



G. The electron effective mass 

In Fig. [8] for selected value of the temperature the 
functions Rc[Z(ti;)] and Im[Z (cj)] are presented {uj S 
(Ojilniax))- On the basis of the plotted courses it can 
be easily noticed, that for the low values of u non-zero is 
only the real part of the wave function renormalization 
factor; this result is directly correlated with the order 
parameter's behavior. For higher frequencies the shape 
of the real and imaginary part of Z iuj) is also plainly 
dependent on the Eliashberg function's form. 

The values of the wave function renormalization factor 
on the complex plane for the selected values of the tem- 
perature is presented in Fig. [HI In particular, in Fig.[S] 
(A) we have shown an exemplary low-temperature solu- 
tion which is characterized by a very complicated course; 
whereas in Fig. [9] (B) we have placed a few selected 
"high" -temperature solutions, which are presenting quite 
similar curves. 

In the framework of the Eliashberg formalism the real 
part of the wave function renormalization factor enables 
the determination of the electron effective mass (?n*). 
In particular, the ratio of to* to the bare electron mass 
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FIG. 8: The dependence of the real and imaginary part of 
the wave function renormalization factor on the frequency for 
T = 17.4 K. The rescaled Eliashberg function is also plotted. 
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FIG. 9: The values of Z {ui) on the complex plane for the 
selected temperatures. The lines with open symbols are for 
Lo £ (0, flmax); the solid lines correspond to uj £ (f2max,i^c). 
(A) Graph for T = 5 K. (B) The courses in the area of the 
"high" temperatures. 



(toe) is given by: TO*/TOe = Re [Z (0)]. Let us notice that 
Re[Z (0)] takes the highest value for T = Tc (similarly 
as Zm=i)- Thus, [to*]jj^3^^ is equal to 2.36TOe. 



III. SUMMARY 

The basic thermodynamic quantities for the supercon- 
ducting phase, that induces in calcium under the influ- 
ence of the pressure 120 GPa, were determined. All cal- 
culations were conducted in the framework of the Eliash- 
berg formalism. In the first step it was shown, that the 
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Coulomb pseudopotential takes very high value equal 
to 0.215. Next it was stated, that the classical Allen- 
Dynes formula significantly underestimates the critical 
temperature. According to the above, the parameters 
in the Allen-Dynes formula were determined once again 
adapting the formula to the strict results predicted by 
the Eliashberg equations. The dimensionless parameter 
Ri = 4.10 was also determined. The obtained result 
means, that the value of the energy gap at the temper- 
ature of zero Kelvin considerably exceeds the value pre- 
dicted by the BCS theory. In the last step, it has been 
stated that the electron effective mass takes the highest 
value equal to 2.36TOe for T = Tc- 

The method of analysis based on the Eliashberg equa- 
tions in the mixed representation allows to determine the 
dependence of the order parameter on the frequency for 
the range from to uic- On the basis of the conducted 
calculations it was stated, that the values of A (w) on the 
complex plane lie on the characteristic spirals with the 



radius decreasing together with the temperature growth. 
The important notification is related to the fact, that 
for u) € (0, Omax) the effective potential of the electron- 
electron interaction is attractive, while for u) € (fimax, Wc) 
is repulsive. 
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